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Abstract 

Optimal Morse matchings reveal essential structures of cell complexes which lead to powerful tools 
to study discrete geometrical objects, in particular discrete 3-manifolds. However, such matchings are 
known to be NP-hard to compute on 3-manifolds, through a reduction to the erasability problem. 

Here, we refine the study of the complexity of problems related to discrete Morse theory in terms 
of parameterized complexity. On the one hand we prove that the erasability problem is iy[P] -complete 
on the natural parameter. On the other hand we propose an algorithm for computing optimal Morse 
matchings on triangulations of 3-manifolds which is fixed-parameter tractable in the treewidth of the 
bipartite graph representing the adjacency of the 1- and 2-simplexes. This algorithm also shows fixed 
parameter tractability for problems such as erasability and maximum alternating cycle-free matching. 
We further show that these results are also true when the treewidth of the dual graph of the triangulated 
3-manifold is bounded. Finally, we investigate the respective treewidths of simplicial and generalized 
triangulations of 3-manifolds. 

Keywords: discrete Morse theory, parameterized complexity, fixed parameter tractability, treewidth, 
completeness, computational topology, collapsibility, alternating cycle-free matching 

1 Introduction 

Classical Morse theory |[37]| relates the topology of a manifold to the critical points of scalar functions 
defined on it, providing efficient tools to understand essential structures on manifolds. Forman ll20l recently 
extended this theory to arbitrary cell complexes. In this discrete version of Morse theory, alternating cycle- 
free matchings in the Hasse diagram of the cell complex, so-called M6>r^^ matchings, play the role of smooth 
functions on the manifold (201 [TTl. For example, similarly to the smooth case ll39l , a closed manifold 
admitting a Morse matching with only two unmatched (critical) elements is a sphere ll20l . The construction 
of specific Morse matchings has proven to be a powerful tool to understand topological ll20ll26l[27l[3n[32l , 
combinatorial EH [25 , 29 1 and geometrical [24, 30, 41 , 40] structures of discrete objects. 

Morse matchings that minimize the number of critical elements are known as optimal matchings lISTIl . 
Together with their number and type of critical elements, these are topological (more precisely simple ho- 
motopy) invariants of the cell complex, just like in the case of the sphere described above. Hence, computing 
optimal matchings can be used as a purely combinatorial technique in computational topology flAj . More- 
over, optimal Morse matchings are useful in practical applications such as volume encoding |[33l|42]l, or 
homology and persistence computation B30l[23l . 

However, constructing optimal matchings is known to be NP-hard on general 2-complexes and on 3- 
manifolds [27l |3T1- This result follows from a reduction to this problem from the closely related 
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erasability problem: how many faces must be deleted from a 2-dimensional simplicial complex before it 
can be completely erased, where in each erasing step only external triangles, i.e. triangles with an edge 
not lying in the boundary of any other triangle of the complex, can be removed 1 19 |? Despite this hardness 
result, large classes of inputs - for which worst case running times suggest the problem is intractable - allow 
the construction of optimal Morse matchings in a reasonable amount of time using simple heuristics |[32l . 
Such behavior suggests that, while the problem is hard to solve for some instances, it might be much easier 
to solve for instances which occur in practice. As a consequence, this motivates us to ask what parameter 
of a problem instance is responsible for the intrinsic hardness of the optimal matching problem. 

In this article, we study the complexity of Morse type problems in terms of parameterized complexity. 
Following Downey and Fellows [15], an NP-complete problem is cdXltA fixed-parameter tractable (FPT) 
with respect to a parameter /c G N, if for every input with parameter less or equal to k, the problem can be 
solved in 0{f{k) • n^^^^) time, where / is an arbitrary function independent of the problem size n. For NP- 
complete but fixed-parameter tractable problems, we can look for classes of inputs for which fast algorithms 
exist, and identify which aspects of the problem make it difficult to solve. Note that the significance of an 
FPT result strongly depends on whether the parameter is (i) small for large classes of interesting problem 
instances and (ii) easy to compute. 

In order to also classify fixed-parameter intractable NP-complete problems, Downey and Fellows fl5\ 
propose a family of complexity classes called the W -hierarchy: FPT C C VF[2] C • • • C VF[P] C 

XP. The base problems in each class of the VF-hierarchy are versions of satisfiability problems with 
increasing logical depth as parameter. Class W[P] contains the satisfiability problems with unbounded 
logical depth. The rightmost complexity class XP of the VF-hierarchy contains all problems which can be 
solved in 0{n^) time where k is the parameter of the problem. 

Here, we use the notion of the VF-hierarchy in a geometric setting. More precisely, we determine the 
hardness of Morse type problems using the mathematically rigid framework of the VF-hierarchy. Our first 
main result shows that the erasability problem is VF[P] -complete (Theorems [s] and |4]), where the parameter 
is the natural parameter - the number of cells that have to be removed. In other words, we prove that the 
erasability problem is fixed-parameter intractable in this parameter. From a discrete Morse theory point 
of view, this reflects the intuition that reaching optimality in Morse matchings requires a global (at least 
topological) context. In this way, we also show that the VF-hierarchy as a purely complexity theoretical tool 
can be used in a very natural way to answer questions in the field of computational topology. Although 
there are many results about the computational complexity of topological problems |[2l[T0l[T9l[35l|43l, to 
the authors' knowledge, erasability is the first purely geometric problem shown to be -complete. 

Our second main result refines the observation that simple heuristics allow us to compute optimal match- 
ings efficiently. For general 2-complexes (and 3-manifolds), the problem reduces directly to finding a maxi- 
mal alternating cycle-free matching on a spine, i.e., a bipartite graph representing the 1- and 2-cell adjacen- 
cies |[3l|27l|33l (Lemma[T]). To solve this problem, we propose an explicit algorithm for computing maximal 
alternating cycle-free matchings which is fixed-parameter tractable in the treewidth of this bipartite graph 
(Theorem [5]). Furthermore, we show that finding optimal Morse matchings on triangulated 3-manifolds is 
also fixed-parameter tractable in the treewidth of the dual graph of the triangulation (Theorem [6]), which is 
a common parameter when working with triangulated 3-manifolds 1 10|. 

Finally, we use the classification of simplicial and generalized triangulations of 3-manifolds to investi- 
gate the "typical" treewidth of the respective graphs for relevant instances of Morse type problems. In this 
way, we give further information on the relevance of the fixed parameter results. The experiments show that 
the average treewidths of the respective graphs of simplicial triangulations of 3-manifolds are particularly 
small in the case of generalized triangulations. Furthermore, experimental data suggest a much more restric- 
tive connection between the treewidth of the dual graph and the spine of triangulated 3-manifolds than the 
one stated in Theorem [6l 
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2 Preliminaries 



Triangulations Throughout this paper we mostly consider simpUcial complexes of dimensions 2 and 3, 
although most of our results hold for more general combinatorial structures. All 2-dimensional simplicial 
complexes we consider are (i) pure, i.e., all maximal simplexes are triangles (2-simplexes) and (ii) strongly 
connected, i.e., each pair of triangles is connected by a path of triangles such that any two consecutive 
triangles are joined by an edge (1-simplex). All 3-dimensional simplicial complexes we consider are tri- 
angulations of closed 3-manifolds, that is, simplicial complexes whose underlying topological space is a 
closed 3-manifold. In particular every 3-manifold can be represented in this way |[36l . We will refer to these 
objects as simplicial triangulations of 3-manifolds. 

In Section [5] we briefly concentrate on a slightly more general notion of a generalized triangulation of a 
3-manifold, which is a collection of tetrahedra all of whose faces are affinely identified or "glued together" 
such that the underlying topological space is a 3-manifold. Generalized triangulations use far fewer tetrahe- 
dra than simplicial complexes, which makes them important in computational 3-manifold topology (where 
many algorithms are exponential time). Every simplicial triangulation is a generalized triangulation, and the 
second barycentric subdivision of a generalized triangulation is a simplicial triangulation ll36l , hence both 
objects are closely related. 

For the remainder of this article, we will often consider 2-dimensional simplicial complexes as part of a 
simplicial triangulation of a 3-manifold. 

Erasability of simplicial complexes Let A be a 2-dimensional simplicial complex. A triangle t E A is 
called external if t has at least one edge which is not in the boundary of any other triangle in A; otherwise 
t is called internal. Given a 2-dimensional simplicial complex A and a triangle t E A, the 2-dimensional 
simplicial complex obtained by removing (or erasing) t from A is denoted by A \ t. In addition, if A' is 
obtained from A by iteratively erasing triangles such that in each step the erased triangle is external in the 
respective complex, we will write A ^ A^ We say that the complex A is erasable if A ^ 0, where in this 
context denotes a complex with no triangle. Finally, for every 2-dimensional simplicial complex A we 
define er(A) to be the size of the smallest subset Aq of triangles of A such that A \ Aq ^ 0. The elements 
of Aq are called critical triangles and hence er( A) is sometimes also referred to as the minimum number of 
critical triangles of A. Determining er(A) is known as the erasability problem 1 19|. 

Problem 1 (Erasability). 

Instance: A 2-dimensional simplicial complex A. 
Parameter: A non-negative integer k. 
Question: Is er(A) < kl 

Hasse diagram and spine Given a simplicial complex A, one defines its Hasse diagram i7 to be a directed 
graph in which the set of nodes of H is the set of simplexes of A, and an arc goes from r to a if and only 
if a is contained in r and dim(cr) + 1 = dim(T). Let Hi C H ho the bipartite subgraph spanned by all 
nodes of H corresponding to i- and i + 1-dimensional simplexes. In particular. Hi describes the adjacency 
between the 2-simplexes and 1-simplexes of A, and will be called the spine of the simplicial complex A. 
The spine of a simplicial complex will be one of the main objects of study in this work. 

Matchings By a matching of a graph G = (A^, A) we mean a subset of arcs M d A such that every node 
of N is contained in at most one arc in M. Arcs in M are called matched arcs and the nodes of the matched 
arcs are called matched nodes. Nodes and arcs which are not matched are referred to as unmatched. By the 
size of a matching M we mean the number of matched arcs. A matching M is called a maximum matching 
of a graph G if there is no matching with a larger size than the size of M. 
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Morse matchings Let H be the Hasse diagram of a simplicial complex A and M be a matching on H. 
Let H{M) be the directed graph obtained from the Hasse diagram by reversing the direction of each arc of 
the matching M. If H{M) is a directed acyclic graph, i.e., H(M) does not contain directed cycles, then M 
is a Morse matching flT]. Furthermore, the number q of unmatched vertices representing z-simplexes of A 
is called the number of critical i-dimensional simplexes and the sum c(M) = ^ - q is said to be the total 
number of critical simplexes. 

The motivation to find optimal Morse matchings is given by the following fundamental theorem of 
discrete Morse theory due to Forman which deals with simple homotopy ||T2|. 

Theorem 1 (|20|). Let M be a Morse matching on a simplicial complex A. Then A is simple homotopy 
equivalent to a CW -complex with exactly one d-cellfor each critical d-simplex of M. 

In other words, a Morse matching with the smallest number of critical simplexes gives us the most 
compact and succinct topological representation up to homotopy. For more information about the basic 
facts of Morse theory we refer the reader to Forman's original work |[20l . This motivates a fundamental 
problem in discrete Morse theory, optimal Morse matching, as a decision problem in the following form. 

Problem 2 (Morse Matching). 

Instance: A simplicial complex A. 
Parameter: A non-negative integer k. 
Question: Is there a Morse matching M with c(M) < kl 

Note that Erasability can be restated as a version of Morse Matching where only the number of 
unmatched 2-simplexes (that is, C2(M)) is counted lISTl . 

Complexity of Morse matchings The compelexity of computing optimal Morse matchings is linear on 
1 -complexes (graphs) [12011 and 2-complexes that are manifolds 1311 . Joswig and Pfetsch [27 J prove that 
if you can solve ERASABILITY in the spine of a 2-simplicial complex in polynomial time, then you can 
solve Morse matching in the entire complex in polynomial time. The proof technique easily extends to 
3-manifolds, leading to the following lemma which has been mentioned in previous works |[33l [3l. 

Lemma 1. Let M be a Morse matching on a triangulated 3-manifold A. Then we can compute a Morse 
matching in polynomial time which has exactly one critical 0-simplex and one critical 3-simplex, such 
that c{M') < c{M). 

In other words, answering ERASABILITY on the spine is the only difficult part when solving MORSE 
Matching on a 3-manifold. In Section |4] we show that if a spine has bounded treewidth, then we can 
solve ERASABILITY in linear time. Lemma [T] therefore generalizes this result to MORSE MATCHING on 
3-manifolds. 

3 W[P]-Completeness of the erasability problem 

In order to prove that ERASABILITY is VF[P] -complete in the natural parameter, we first have to take a closer 
look at what has to be considered when proving hardness results with respect to a particular parameter. 

Definition 1 (Parameterized reduction). A parameterized problem L reduces to a parameterized problem 
V , denoted by L <fpt L', if we can transform an instance (x, k) of L into an instance 9{k)) of V in 
time (where f and g are arbitrary functions), such that (x, k) is a yes-instance of L if and only 

if{x\ g{k)) is a yes-instance of V . 
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As an example, Egecioglu and Gonzalez Ull reduce Set Cover to Erasability to show that Erasabil- 
ITY is NP-complete. Since their reduction approach turns out to be a parameterized reduction, these results 
can be restated in the language of parameterized complexity as follows. 

Corollary 1. Set Cover <fpt Erasability, therefore Erasability is W[2]-hard. 

This shows that, if the parameter k is simultaneously bounded in both problems, ERASABILITY is at 
least as hard as Set Cover. In this section we will determine exactly how much harder Erasability is 
than Set Cover, which is [2] -complete. Namely, we will show that Erasability is VI/[P] -complete 
in the natural parameter k. This will be done by i) using a VK[P] -complete problem as an oracle to solve 
an arbitrary instance of ERASABILITY (Theorem [sj which shows that ERASABILITY is in VF[P]), and ii) 
reducing an arbitrary instance of a suitable problem which is known to be VF[P] -complete to an instance of 
Erasability (Theoremji} which shows that Erasability is VI/[P]-hard). 

There are only a few problems described in the literature which are known to be VF[P] -complete |[T6l 
p. 473]. Amongst these problems, the following is suitable for our purposes. 

Problem 3 (Minimum Axiom Set). 

Instance: A finite set S of sentences, and an implication relation R consisting 

of pairs ((7, s) where U C S and s ^ S. 
Parameter: A positive integer k. 

Question: Is there a set Sq C S (called an axiom set) with l^ol < and a 
positive integer n, for which Sn = S, where we define Si,l < i < 
n, to consist of exactly those 5 E 5 for which either s G Si-i or 
there exists a set ?7 C Si-i such that (?7, s) G P? 

Theorem 2 (CSl). Minimum Axiom Set is W[P]-complete. 

In this paper, we show that, preserving the natural parameter k. Minimum Axiom Set is both at least 
and at most as hard as Erasability. 

Theorem 3. Erasability <fpt Minimum Axiom Set, therefore Erasability is in W[P]. 

Theorem [s] shows that ERASABILITY is at most as hard as the hardest problems in VF[P]. Please refer 
to the full version of this paper for a detailed proof of Theorem |3] 

In order to show that it is in fact amongst the hardest problems in this class we first need to build some 
gadgets. 

Definition 2 (Gadgets for the hardness proof of Erasability). Let {S, P, k) be an instance of Minimum 
Axiom Set. 

Let s ^ S be a sentence. By an 5-gadget or sentence gadget we mean a triangulated 2-dimensional 
sphere with 2n-\- m punctures as shown in Figure^ where m is the number of relations ((7, s) ^ R and n 
is the number of relations (U^u) ^ R such that s ^ U. 

Let ((7, s) ^ R be a relation. A (U, 5)-gadget or implication gadget is a collection of\U\ + 1 sentence 
gadgets for each sentence ofUU{s} together with 2\U\ nested tubes as shown in Figure^such that (i) two 
tubes are attached to two punctures of the u-gadgetfor each u G U and (ii) all2\U\ boundary components 
at the other side of the tubes are identified at a single puncture of the s-gadget. 

Then, by construction the following holds for the (?7, 5)-gadget. 

Lemma 2. A (U^ s)- gadget can be erased if and only if all sentence gadgets corresponding to sentences in 
U are already erased. 
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Figure 1: Example of a sentence gadget with m = 2 relations (U^s) and n = 3 relations (U^u) with 
additional tubes. 




Proof. Clearly, if all sentence gadgets corresponding to sentences in U are erased, the whole gadget can be 
erased tube by tube. If, on the other hand, one of the sentence gadgets still exists, this gadget together with 
the two tubes connected to it build a complex without external triangles which thus cannot be erased. □ 

With these tools in mind we can now prove the main theorem of this section. 

Theorem 4. Minimum Axiom Set <fpt Eras ability, even when the instance of Eras ability 
is a strongly connected pure 2 -dimensional simplicial complex A which is embeddable in R^. Therefore 
Eras ABILITY is W[P]-hard. 

The simplicial complex A (Figure [s]) constructed to prove -hardness of ERAS ABILITY is in fact 
embeddable into R^. This means that, even in the relatively well behaved class of embeddable 2-dimensional 
simplicial complexes, ERASABILITY when bounding the number of critical simplexes is still likely to be 
inherently difficult. Please refer to the full version of this paper for a detailed proof of Theorem |4j The 
VF[P] -completeness result implies that if ERASABILITY turns out to be fixed parameter tractable, then 
W[P] = FPT, i.e., every problem in W[P] including the ones lower in the hierarchy would turn out to 
be fixed parameter tractable, an unlikely and unexpected collapse in parameterized complexity. Also, it 
would imply that the n- variable SAT problem can be solved in time 2^^^^ , that is, better than in a brute force 
search HI. With respect to this result, if we want to prove fixed parameter tractability of ERASABILITY, the 
parameter must be different from the natural parameter. 



4 Fixed parameter tractability in the treewidth 

In this section, we prove that there is still hope to find an efficient algorithm to solve MORSE MATCHING. 
We give positive results for the field of discrete Morse theory by proving that ERASABILITY and MORSE 
Matching are fixed parameter tractable in the treewidth of the spine of the input simplicial complex, and 
also in the dual graph of the problem instance in case it is a simplicial triangulation of a 3-manifold. 
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Figure 3: Constructing an instance of ERASABILITY from an instance (5, R) of Minimum Axiom Set, 
where S = {a, 6, c, d, e, /, g, h, i} and i? = {({c, rf, e}, i), ({/, g, h}, z), ({6}, c), ({a, rf}, ^)}. 




Figure 4: Example of a nice tree decomposition (left) of the spine of a non-manifold 2-dimensional simpli- 
cial complex (right). 

4.1 Treewidth 

Definition 3 (Treewidth). A tree decomposition of a graph G is a tree T together with a collection of bags 
{Xi], where i is a node ofT. Each bag Xi is a subset of nodes ofG, and we require that (i) every node of 
G is contained in at least one bag Xi (node coverage J; (ii)for each arc ofG, some bag Xi contains both its 
endpoints fare coverage); and for all bags Xi, Xj and Xj^ ofT, if Xj lies on the unique simple path from 
Xi to Xk in T, then XiH Xj. ^ Xj (coherence! 

The width of a tree decomposition is defined as max \ Xi\ — 1, and the treewidth ofG is the minimum 
width over all tree decompositions. We will denote the treewidth ofG by tw(G). 

For bounded treewidth, computing a tree decomposition of a graph G = (y, E) of width < k has running 
time 0{f{k)\V\) [5| due to an algorithm by Bodlaender. Regarding the size of f{k)'. using the improved 
algorithm by Perkovic and Reed |38 1, at most 0(/c^) recursive calls of Bodlaender's improved linear time 
fixed-parameter tractable algorithm for bounded treewidth from f6l are needed. This latter algorithm in turn 
is said to have a constant factor f{k) which is "at most singly exponential in fc". For further reading on the 
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running times of tree decomposition algorithms see ll4l[28ll. 

Definition 4 (Nice tree decomposition). A tree decomposition [Xi \ i ^ I^T) is called a nice tree decompo- 
sition if the following conditions are satisfied: 

1. There is a fixed bag with | = 1 acting as the root ofT( in this case is called the root bag J. 

2. If bag Xj has no children, then \Xj \ — 1 (in this case Xj is called a leaf bag| 

3. Every bag of the tree T has at most two children. 

4. If a bag Xi has two children Xj and X^, then Xi — Xj — X^ (in this case Xi is called a join bag J. 

5. If a bag i has one child j, then either 

(a) = + 1 and Xj C Xj (in this case Xi is called an introduce bag J, or 

(b) |Xj I = + 1 and Xi C Xj (in this case Xi is called a forget bagj. 

A given tree decomposition can be transformed into a nice tree decomposition (Figure |4]) in linear time: 

Lemma 3 (|28l). Given a tree decomposition of a graph G of width w and 0{n) bags, where n is the 
number of nodes ofG, we can find a nice tree decomposition ofG that also has width w and 0{n) bags in 
time 0{n). 

4.2 Alternating cycle-free matchings 

Given a graph G = (A^, A) and a matching M C A on G, an alternating path is a sequence of pairwise 
adjacent arcs such that each matched arc in the sequence is followed by an unmatched arc and conversely 
An alternating cycle of M is a closed alternating path. Matchings which do not have any such alternating 
cycle are called alternating cycle-free matchings. From the definition of Morse matching, we can state 
Erasability in the language of alternating cycle-free matchings as follows: 

Problem 4 (Alternating cycle-free matching). 
Instance: A bipartite graph G = (A^i U A^2, 
Parameter: A nonnegative integer k. 

Question: Does G has an alternating cycle-free matching M with at most k 
unmatched nodes in A^i? 

Specifically, if G = i^i is the spine for some simplicial complex A, then ERASABILITY is equivalent to 
the Alternating cycle-free matching problem. 

4.3 FPT algorithm for the alternating cycle-free matching problem 

Courcelle's theorem L13J can be used to show that Alternating cycle-free matching is fixed param- 
eter tractable (please refer to the full version of this paper). However, this is a purely theoretical result, since 
the stated complexity contains towers of exponents in the parameter function. This is the reason why, for the 
remainder of this section, we focus on the construction of a linear time algorithm to solve ALTERNATING 
CYCLE-FREE MATCHING for inputs of bounded treewidth with a significantly faster running time. 

Theorem 5. Let G — {Ni U N2,A) be a simple bipartite graph with a given nice tree decomposition 
{Xi I i G /, T). Then the size of a maximum alternating cycle-free matching of G can be computed in 
0(4^ ' ' \og{w) • n) time, where n — \N\ and w denotes the width of the tree decomposition. 
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Algorithm overview Our algorithm constructs alternating cycle-free matchings of G along the nice tree 
decomposition [Xi | z G /, T) of G, from the leaves up to the root, visiting each bag exactly once. In the 
following we will denote by F^, the set of nodes which are already processed and forgotten by the time Xi 
is reached; we call this the set of forgotten nodes. At each bag Xi of the decomposition, we construct a set 
M{i) representing all valid alternating cycle-free matchings in the graph induced by the nodes in U F^. 

The leaf bags contain a single node of G, and the only matching is thus empty. At each introduce bag 
Xi = Xj U {x}, each matching M of M{j) can be extended to several matchings as follows. The newly 
introduced node x can be either left unmatched, or matched with one of its neighbors as long as it generates 
a valid and cycle-free matching with M. At each join bag Xi = Xj = X/^, M{i) is build from the valid 
combinations of pairs of matchings from M ( j) and M{k). The final list of valid matchings is then evaluated 
at the root bag r. 

However, this final list M{r) contains an exponential number of matchings. Fortunately, the nice tree 
decomposition allows us to group together, at each step, all matchings M that coincide on the nodes of X^. 
Indeed, the algorithm takes the same decisions for all the matchings of the group. We can thus store and 
process a much smaller list M{j) of matchings containing only one representative M of each group. In 
each group, we choose one with the smallest number of unmatched nodes so far. This grouping takes place 
at the forget and join bags. This makes the algorithm exponential in the bag size, not the input size. The 
algorithm is described step-by-step and illustrated in the full version of this paper. 

Matching data structure The structure storing an alternating cycle-free matching M in a set M{i) must 
be suitable for checking the matching validity whenever a matching is extended at an introduce bag or a join 
bag. It must store which nodes are already matched in M to avoid matching a node of G twice {matching 
condition). We use a binary vector v(M), where the x-th coordinate is 1 if node x G X^ is matched and 
otherwise. Checking the matching condition and updating when nodes are matched has thus a constant 
execution time 0(1). 

Also, the structure must store which nodes are connected by an alternating path in M to avoid closing 
a cycle when extending or combining M {cycle-free condition). When matching two nodes x and y, an 
alternating cycle is created if there exists an alternating path from a neighbor of x to a neighbor ofy. To test 
this, we use a union-find structure [44J uf (M), storing for each matched node x the index of a matched node 
c{x) connected to x by an alternating path in M. For a subset of matched nodes which are all connected to 
each other, the component index c is chosen to be the node with the lowest index. For each unmatched node, 
we store the ordered list of component indexes of neighbor matched nodes. The cycle-free condition check 
reduces to find calls on the adjacent lists, and the update of the structure when increasing the matching size 
reduces to union calls, both executing in near-constant time. All the matchings are stored in a hash structure 
to allow faster search for duplicates. Finally, we can return not only the maximal cycle-free matching size, 
but an actual maximal cycle-free matching by storing, along with each representative matching, a binary 
vector of size |Xi U F^l with all the matched nodes so far. 

Grouping Traversing the nice tree decomposition in a bottom-up fashion, each node appears in a set of 
bags that form a subtree of the tree decomposition {coherence requirement). This means that, whenever a 
node is forgotten, it is never introduced again in the bottom-up traversal. 

A naive version of the algorithm described above would build the complete list of valid alternating 
cycle-free matchings: the set M{i) would contain all valid matchings in the graph induced by the nodes in 
Xi[JFi. In particular, for each matching M ^ M{i) the algorithm would store the binary vector v(M) and 
the union-find structure uf (M) on Xi U Fi. However, it is sufficient to store the essential information about 
each M by restricting the union-find structure uf (M) and the binary vector v(M) only to the nodes in the 
bag Xi (for any matched node x G X^, node c{x) of the union-find structure is then chosen inside Xi). 



9 



n 


# triangulations 


min w 


max w 


w 


minw 


max w 


w (dual) 


1 


4 


(3) 


1 


2 


1.50 


(1.67) 








0.00 


2 


17 


(12) 


1 


3 


2.47 


(2.42) 


1 


1 


1.00 


3 


81 


(63) 


1 


3 


2.51 


(2.49) 


1 


2 


1.60 (1.52) 


4 


577 


(433) 


1 


5(4) 


2.77 


(2.73) 


1 


3 


1.91 (1.87) 


5 


5184 


(3961) 


1 


6(5) 


2.95 




1 


4 


2.16 (2.18) 


6 


57753 


(43584) 


1 


6 


3.16 


(3.19) 


1 


4 


2.31 (2.35) 


7 


722765 


(538409) 


1 


7 


3.35 


(3.40) 


1 


4 


2.45 (2.50) 



Table 1 : Treewidths of the spine (left) and of the dual graphs (right) of closed generalized triangulations up 
to 7 tetrahedra. The values in brackets are for 1- vertex triangulations. 



More precisely, we define an equivalence relation on the matchings of such that M M' if and 
only if v(M) = v(M') and uf (M) = uf (M^) on the nodes of Xi. Since two equivalent matchings only 
differ on the forgotten nodes F^, the validation of the matching and cycle-free conditions of any extension 
of M or M' (or any combination with a third equivalent matching M") will be equal from now on. 

Since we are interested in the alternating cycle-free matching with the minimum number of unmatched 
nodes, for each equivalence class we will choose a matching M with the minimum number m{M) of 
unmatched forgotten nodes as class representative. This number m{M) is stored together with (v,uf) 
for each equivalence class of M{i) = M{i)/ ^i. In addition, we can compute the alternating cycle-free 
matching of maximum size by storing the complete binary vector v along with m{M) (since the matching 
is cycle-free, this is sufficient to recover the set of arcs defining the matching). 



Execution time complexity To measure the running time we need to bound the number of equivalence 
classes of M{i). Let wi be the number of nodes in Xi. The number of equivalence classes of M{i) is then 
bounded above by the number of possible pairs (v, uf ) on wi nodes. The union-find stores for each node 
X, either a component node c{x) G or a list of at most wi component nodes, leading to at worst 2^^ 

2 

different lists per node, giving 2^^ possible combinations of lists. Also there are 2^^ possible binary vectors 
V of length wi, therefore there are at worst 2^^ 2^^ elements in (this enumeration includes invalid 

matchings and incoherences between v and uf ). 

The time complexity is dominated by the execution at the join bag where pairs of equivalences classes 
from M{j) and M{k) have to be combined. Therefore we must square the number of equivalence classes 
in each set: the complexity for a join bag is 0(4^ - - \og{w)) (please refer to the full version of this 
paper for details). Since there are 0{n) bags in a nice tree deco mpo sition, the total execution time is in 



0(4^ +w; . ^3 . iog(it;) • n). Finally, as already stated in Section 4.1 for bounded treewidth computing a 
tree decomposition and a nice tree decomposition is linear. So the whole process from the bipartite graph 
to the resulting maximal alternating cycle-free matching is fixed-parameter tractable in the treewidth. Note 
that neither the decomposition nor the algorithm use the fact that the graph is bipartite. 



4.4 Correctness of the Algorithm 

We must check that the algorithm, without the grouping, considers every possible alternating cycle-free 
matching in G and that the grouping occurring at the forget and join bags does not discard the maximal 
matching. 

The node coverage and arc coverage properties of nice tree decompositions (Definition [3]) ensure that 
each node is processed and each arc is considered for inclusion in the matching at one introduce node. Since 
the introduce node discards only matchings that violate either the matching or the cycle condition, and these 
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n 


^ triangulations 


min w 


max w 


w 


min w 


max w 


w (dual) 


5 


1 


6 


6 




6.00 


4 


4 


4.00 


6 


2 


< 7 


< 8 


< 


7.50 


4 


5 


4.50 


7 


5 


< 8 


< 11 


< 


9.40 


4 


6 


5.00 


8 


39 


< 8 


< 14 


< 


11.23 


4 


7 


5.74 


9 


1297 


< 8 


< 18 


< 


13.55 


4 


9 


7.01 


10 


249015 


< 8 


< 22 


< 


16.33 


< 4 


< 13 


< 8.99 



Table 2: Upper bounds and exact values for the treewidths of the spine (left) and of the dual graph (right) of 
simplicial triangulations of 3-manifolds up to 10 vertices. 

violations cannot be legalized by further extensions or combinations of the matchings, all possible valid 
matchings are considered. 

Now, consider two matchings M and M' that are grouped together and represented by M at a forget 
or join bag Xi. In the further course of the algorithm, the representative M is then extended or combined 
with other matchings to form new valid matchings M^ The coherence property of Definition [s] assures 
that no neighbor of a newly introduced node can be a forgotten node, so the extension or combination only 
modifies matchings M and M on nodes of X^, which are represented in the structure of M' . Hence, the valid 
matchings M' actually represent all the valid extensions and combinations of M and M. The grouping thus 
generates all valid and relevant representatives of matchings in order to find a maximal alternating cycle- 
free matching. Moreover, in case M and M' are equivalent and both with the lowest number of forgotten 
unmatched nodes, choosing M or M' as representative leads to the exact same extensions and combinations. 

Finally, let be the alternating cycle-free matching of maximum size of G. In each bag the corre- 
sponding matching must be one of the matchings with the lowest number of unmatched nodes within its 
equivalence class G M{i). Otherwise, a matching in the same class M^, extended and combined as 
Mjn in the sequel of the algorithm would give rise to a matching with fewer unmatched nodes. Therefore, 
the choice of the representative at the forget and join bags never discards the future alternating cycle-free 
matching of maximum size. 

4.5 Treewidth of the dual graph 

Up to this point, we have been dealing primarily with simplicial complexes and their spines. We now turn 
our attention to simplicial triangulations of 3-manifolds and a more natural parameter associated to them. 

Definition 5 (Dual graph). The dual graph of a simplicial triangulation of a 3-manifold T, denoted r(T), 
is the graph whose nodes represent tetrahedra of T, and whose arcs represent pairs of tetrahedron faces 
that are joined together 

We show that, if the treewidth of the dual graph is bounded, so is the treewidth of the spine, as stated by 
the following theorem (please refer to the full version of this paper for the proof). 

Theorem 6. Let G be the spine of a simplicial triangulation of a 3-manifold T. Iftv7{V(T)) < k, then 
tv^(G) < 10k + 9. 

5 Experimental Results 

In Section |3] we have seen that the problem of finding optimal Morse matchings is hard to solve in general. 
In Section]?] on the other hand we proved that in the case of a small treewidth of the spine of a 2-dimensional 
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complex or, equivalently, in the case of a bounded treewidth of the dual graph of a simplicial triangulation of 
a 3-manifold, finding an optimal Morse matching becomes easier. Up to a certain scaling factor, the results 
stated in Section |4] hold for generalized triangulations as well (also, note that the notion of a spine or the 
dual graph can be extended in a straightforward way to generalized triangulations). 

Given this situation, a natural question to ask is the following: What is a typical value for the treewidth 
of the respective graphs of (i) small generic generalized triangulations of 3-manifolds, and (ii) small generic 
simplicial triangulations of 3-manifolds? 

In a series of computer experiments we computed the treewidth of the relevant graphs (i.e., the spine 
and the dual graph) of all closed generalized triangulations of 3-manifolds up to 7 tetrahedra |8|, and all 
simplicial triangulations of 3-manifolds up to 10 vertices lf34l . The computer experiments were done using 
LibTW [45j to compute the treewidth / upper bounds for the treewidth, with the help of the GAP package 
simpconnp |[T71[T8l and the 3-manifold software Regina ||7l|9l. We report the minimal, maximal and 
average treewidths of all triangulations with the same number of vertices in Tables [T] and [2j 

Regarding the treewidth of generalized triangulations of 3-manifolds, we observe that there is a large 
difference between the average treewidth and the maximal treewidth for both the dual graph and the spine. 
In particular, the average treewidth appears to be relatively small. Moreover, there is only a slight difference 
between the data for general closed triangulations and 1-vertex triangulations. This fact is somehow in 
accordance with our intuition since the number of 0-dimensional simplexes should neither directly affect 
the spine nor the dual graph of a generalized triangulation. 

On the other hand, the gap between the maximum treewidth and the average treewidth in the case of sim- 
plicial triangulations of 3-manifolds is relatively small compared to the data for generalized triangulations. 
In addition, the treewidth of the spines of some particularly interesting 2-dimensional simplicial complexes 
(reported in the full version of this paper) is significantly smaller than the (upper bound of the) treewidth of 
simplicial triangulations of 3-manifolds. At this point it is important to note that, while the data concerning 
the spines for simplicial complexes only consists of upper bounds, experiments applying the algorithm for 
the upper bound to smaller graphs and then computing their real treewidths suggest that these upper bounds 
(in average) are reasonably close to the exact treewidth. 

Further analysis shows that the average treewidth of the spines for both generalized and simplicial tri- 
angulations of 3-manifolds is mostly less than twice the treewidth of the dual graph, and hence much below 
the theoretical upper bound given by Theorem [6j Also, the ratio between these two numbers appears to be 
more or less stable for all values shown in Tables [T] and [2j This can be seen as experimental evidence that 
for triangulated 3-manifolds the treewidth of the dual graph is responsible for the inherent difficulty to solve 
Erasability and related problems. 

Despite the small values of n in our tables, there are theoretical reasons to believe that the patterns of 
small treewidth will continue for larger n. For instance, the conjectured minimal triangulations of Seifert 
fibered spaces over the sphere have dual graphs with 0(1) treewidth for arbitrary n. Moreover, following 
recent results of Gabai et al. [21] there are reasons to believe that large infinite classes of topological 
3-manifolds admit triangulations whose treewidths are below provable upper bounds. Investigating these 
upper bounds is work in progress. 
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Appendix 



A Proof of Lemma [T] 

The proof of Lemma [T] actually follows directly from Joswig and Pfetsch's proof of the following lemma. 

Lemma 4 ( ||27]| ). Let M be a Morse matching on a simplicial complex A. Then we can compute a Morse 
matching M' in polynomial time which has exactly one critical 0-simplex, the same number of critical 
simplexes of dimension greater or equal 2 as M, and c{M') < c(M). 

The proof builds a Morse matching from a spanning tree of the primal graph, i.e. the graph obtained 
considering only the vertices and edges of A. For a 3-manifold A, the proof of the previous lemma can be 
applied exactly the same way on the dual graph of A to obtain the following result. 

Lemma 5. Let M be a Morse matching on a closed triangulated 3-manifold A. Then we can compute a 
Morse matching M' in polynomial time which has exactly one critical 3-simplex, the same number of critical 
simplexes of dimension less or equal 1, and c{M') < c(M). 

Since the proof works independently on the primal and dual graph, Lemma [T] is a combination of these 
results. Here, we simply reproduce the proof of Joswig and Pfetsch L27.I verbatim applying it to 3-manifold 
complexes, using Poincare's duality. 

First consider a Morse matching M for a connected 3-manifold A. Let 7(M) be the graph obtained 
from the primal graph of A by removing all arcs (edges of A) matched with triangles in M and let 7* (M) 
be obtained from the dual graph of A by removing all the arcs (triangles in A) where the corresponding 
triangles are matched with edges of A in M. Note that 7(M) and 7* (M) contain all vertices and tetrahedra 
of A, respectively. 

Lemma 6. The graph 7(M) and dual graph 7* (M) are connected. 

Proof Suppose that 7(M) is disconnected. Let N be the set of nodes in a connected component of 7(M), 
and let C be the set of cut edges, that is, edges of A with one vertex in N and one vertex in its complement. 
Since A is connected, C is not empty. By definition of 7(M), each edge in C is matched to a unique 
2-simplex. 

Consider the directed subgraph D of the Hasse diagram consisting of the edges in C and their matching 
2-simplexes. The standard direction of arcs in the Hasse diagram (from the higher to the lower dimensional 
simplexes) is reversed for each matching pair of M, i.e., is a subgraph of H{M). We construct a directed 
path in D as follows. Start with any node of D corresponding to a cut edge ei. Go to the node of D 
determined by the unique 2-simplex ti to which ei is matched to. Then ti contains at least one other cut 
edge 62, otherwise ei cannot be a cut edge. Now iteratively go to 62, then to its unique matching 2-simplex 
T2, choose another cut edge 63, and so on. We observe that we obtain a directed path ei, ri, 62, T2, • • • in 
D, i.e., the arcs are directed in the correct direction. Since we have a finite graph at some point the path 
must arrive at a node of D which we have visited already. Hence, D (and therefore also H{M)) contains a 
directed cycle, which is a contradiction since M is a Morse matching. 

To prove that 7* (M) is connected, we repeat the proof above on the dual graph. □ 

Lemma 1 ( |[33l[27l [3l). Let M be a Morse matching on a triangulated 3-manifold A. Then we can compute a 
Morse matching M' in polynomial time which has exactly one critical 0-simplex and one critical 3-simplex, 
such that c{M') < c{M). 
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Proof. Since 7(M) and 7*(M) are connected, they both have spanning trees, and we will use them to build 
the Morse matching. First pick an arbitrary node ri and any spanning tree of 7(M) and direct all arcs 
away from ri. Then pick an arbitrary tetrahedron (a node in the dual graph) r2 and any spanning tree of 
7*(M) and direct all triangles (arcs in dual graph) away from r2. This yields a maximum Morse matching 
on 7(M) and 7*(M). It is easy to see that replacing the part of M on 7(M) and 7*(M) with this matching 
yields a Morse matching. This Morse matching has only one critical vertex (the root ri) and one critical 
tetrahedron (the root r2). Note that every Morse matching in a triangulated 3-manifold contains at least one 
critical vertex and at least one critical tetrahedron; this can be seen from Theorem [T| Furthermore, the total 
number of critical simplexes can only decrease, since we computed an optimal Morse matching on 7(M) 
and7*(M). □ 
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B Proof of Theorems] 



Theorem 3. Erasability <fpt Minimum Axiom Set, therefore Erasability is in W[P]. 

Proof. We show membership of ERASABILITY in W[P] by reducing a given instance (A, k) of ERASABIL- 
ITY to an instance (5, R, k) of Minimum Axiom Set. 

W. 1. o. g. we can assume that the 2-dimensional simplicial complex A has no external edges (if A has 
external edges we first remove these edges until no external edge exists and reduce the remaining problem 
instance to an instance of Minimum Axiom Set). We now identify the set of triangles of A with the set 
of sentences 5 in a one-to-one correspondence. For every edge e G A we denote the set of all triangles 
containing e by star a (e) C A, we write for the corresponding set of sentences Se C S, and we define the 
set of implication relations R by the relations 

{Se\{s},s) 

for each triangle s ^ Se for all edges e E A. Note that A has no external edges and thus \ {5} 7^ for 
all e. 

In a next step, we show that for all axiom sets Sq C S of size k we have A \ Aq ^ for the associated 
subset of triangles Aq C A of size k. To see that this is true, note that for the augmenting sequence 
So C Si C ... C Sn = S of S, their corresponding subsets of triangles Aq C Ai c . . . C A^ = A and 
i G {1, . . . , n} fixed, all sentences s ^ Si \ Si-i have to occur in a relation {Se \ {5}, s) for some edge e 
with Se \ {s} C Si-i. For the triangle t G A corresponding to s this means that, star a (e) \ t C A^-i. 
Thus, if we assume that all triangles in A^_i are already erased, t must be external and thus can be erased as 
well. The statement now follows by the fact that for z = 1, all triangles in Aq are already erased in A \ Aq 
and hence A \ Aq 0. 

Conversely, let Aq C A be of size k such that A \ Aq ^ 0. Since A has no external triangles but 
A \ Aq ^ 0, there must be external triangles t G A \ Aq and hence for 5 G 5 being the sentence 
corresponding to the triangle t there is a relation (Se \ {5}, s) with Se \ {s} C Sq, where Sq is the set 
of sentences corresponding to the set of triangles Aq. We then define Si to be the union of Sq with all 
sentences s of the type described above and iterating this step results in a sequence of subsets Sq C Si C 
... C Sn = S for some n what proves the result. □ 
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C Proof of Theorem HI 



Theorem 4. Minimum Axiom Set <fpt Eras ability, even when the instance of Eras ability 
is a strongly connected pure 2 -dimensional simplicial complex A which is embeddable in R^. Therefore 
Eras ABILITY is W[P]-hard. 

Proof To show VF[P] -hardness of ERAS ABILITY, we will reduce an arbitrary instance (5, i?, k) from Min- 
imum Axiom Set to an instance (A, k) of Erasability. In order to do so, we will use a sentence gadget 
for each element of S and an implication gadget for each relation R (cf. Definition [2]) to construct a 2- 
dimensional simplicial complex A with a polynomial number of triangles in the input size. 

By construction, we can glue all sentence and implication gadgets together in order to obtain a simplicial 
complex A without any exterior triangles. Note that the only place where A is not a surface is at the former 
m boundary components per sentence gadget corresponding to the m relations in R with the respective right 
hand side. 

For any axiom set Sq C S of size k, let Aq be a set of k triangles, one from each sentence gadget 
corresponding to a sentence in Sq. It follows by Lemma [2| that A \ Aq can be erased to a complex where 
all the sentence gadgets s corresponding to relations {U,s),U C Sq, have external triangles. Since Sq is an 
axiom set, iterating this process erases the whole complex A. 

Conversely, let Aq be a set of k triangles such that A \ Aq 0. First, note that erasing a triangle of any 
tube of an implication gadget always allows us to remove the sentence gadget at the right end of this tube. 
Hence, w. 1. o. g. we can assume that all k triangles in Aq are triangles of some sentence gadget in A. Now, 
if any sentence gadget contains more than one triangle of Aq we delete all additional triangles obtaining a 
set Aq of k' triangles, k' < k, such that A \ Aq --^ and thus the corresponding set of sentences is an axiom 
set of size k' <k. 

The result now follows by the observation that A can be realized by at most a quadratic number of 
triangles in the input size of (5, R^k). □ 
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D Fixed parameter tractability of Alternating cycle-free matching 
from Courcelle's theorem 



Courcelle's celebrated theorem fTSl states that all graph properties that can be defined in Monadic Second- 
Order Logic (MSOL) can be decided in linear time when the graph has bounded treewidth. Here, we want 
to use Courcelle's theorem to show that problems in discrete Morse theory are fixed parameter tractable 
in the treewidth of some graph associated to the problem. However, it is not obvious how to directly state 
Erasability and Morse Matching in MSOL. Instead, we will apply Courcelle's theorem to Alter- 
nating CYCLE-FREE MATCHING which by the comment made in Section|42]is a graph theoretical problem 
equivalent to ERASABILITY. 

Theorem 5. Let w >l. Given a bipartite graph with tw{G) < w, ALTERNATING CYCLE-FREE MATCH- 
ING can be solved in linear time. 

Proof. We will write a MSOL formulation of ALTERNATING CYCLE-FREE MATCHING based on the fact 
that M is an alternating cycle-free matching if and only if M is a matching and every induced M-subgraph 
contains a node of degree 1 Ii22il : 

maxM : Vx G A^[-i3ai, a2 G M(ai 7^ a2 A inc(x^ ai) A inc(x^ <^2))] 
AVM' C M(3a G M',3x G N[inc(x, a)A(\fxi G 7V(^3ai G M\x ^ xAadj{x,xi)Ainc{xi,ai))))]) 

where inc{x^ a) is the incidence predicate between node x and arc a and adj{x^ x') is the adjacency predicate 
between node x and node x' . The above statement can be translated to plain English as follows: "Find the 
largest matching M of G, where each node is incident to at most one arc, such that in every subset M' of 
the matching M there exists a matched node x in M' such that its only neighbor matched in M' is the other 
endpoint of the unique matched arc incident to x. □ 
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E Algorithm for Alternating cycle-free matching: step by step 



The algorithm visits the bags of the nice tree decomposition bottom-up from the leaves to the root evaluating 
the corresponding mappings in each step according to the following rules (Figure [8]). 

Leaf bag The set of matchings M{i) of a leaf bag Xi = {x} is trivial with a unique empty matching M 
represented by v(M) = [0], and uf(M)(x) defined by as an empty list, associated with m{M) = 0. 

Introduce bag Let = Xj U {x} be an introduce bag with child bag Xj. The set of valid matchings 
M{i) is built from M{j) by introducing x in each matching M G M{j), generating several possible 
matchings M^ We can always introduce x as an unmatched node, then M is extended on x by setting 
v{M')\x = and updating uf(M') with the ordered list of components for each matched neighbor of x. 
In addition, for each unmatched neighbor y e Xj, we can introduce x as a matched node in the following 
way. We match both x and y in M and set ^{M')^^ = 1 and v(M')|^ = 1. If the intersection of the list of 
neighbor components of x and y is empty, then the matching of x and y does not create a cycle. In this case 
is a valid extension of M. The update of the union-find structure must then reflect the extensions of 
all alternating paths through arc {x, y}. We perform in uf(M') a union operation for x and all its matched 
neighbors (including y), and for y and all its matched neighbors. We also add the merged component index 
c{x) to the list of neighbor components of each unmatched neighbor of x and y. Then we include all valid 
extensions to M{i), reducing v(M') by calling find for each node and neighbor component list entry, 
and we set m{M) — m{M') for all extensions M' of M. 




Figure 5: Details of the decisions at an introduce bag. 



Running time There are at most 2^^ extended matchings M' for bag Xi (including all invalid ones), 
where it;^ = = |Xj | + 1 (a new possible matching can be generated only once). Each new matching is 
validated by a direct lookup at v(M') and ordered list comparison, leading to a linear time Wi. The update 
of each structure requires constant time for each matched neighbor of x and almost linear time 0{wi) plus 
the sorted insertion 0{wi • \og{wi)) for each unmatched neighbor, and there are at most Wi neighbors in the 
bag. Thus, the total running time of an introduce bag is in 0(2^* • wf • log{wi)). 

Forget bag Let Xi = Xj \ {x} be a forget bag with child bag Xj 3 x. While the set of all possible 
matchings on Xi U Fi does not change (M{j) = M{i)), the equivalence relation possibly identifies 
more matchings than ^j. For each matching M G yW(j), a new matching is obtained by deleting 
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coordinate x of v(M). If c(x) = x, uf(M) needs to be updated. To do so, the set of nodes Xi is traversed 
twice, once to look for node y ^ x of minimal index such that c{y) = c{x) (eventually, y is empty), and 
a second time to replace xhy y each time x is used as a component index. If x was unmatched in M (i.e., 
v(M)|^ 0), then we set m{M^) m{M) + 1, otherwise we set m{M') m{M). Once the SQiM{j)' of 
all the newly generated M' is computed, M{i) is obtained as the quotient of M{jy by the equivalence 
relation on Xi. More precisely, each pair (M^, M") G M{jy'^ is tested for equality on both v and uf. If 
they are equal, the one with the lowest m is defined to be the new representative in M{i). 
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forget unmatched node: ++m 



forget unmatched node: ++m 




Figure 6: Details of the decisions at a forget bag. 



Running time Each new matching M' is obtained from a single element of in worst-case time 

0{wl • \og{wi)). Equivalent matchings are detected on-the-fly when filling the hash structure of M{i), and 
each equivalence test is linear in wf. The complexity is thus in 0(2^^^^^ • w'j • log{wj)). 




Figure 7: Details of the decisions at a join bag. 



Join bag Let Xi = Xj = X/^ be a join bag with child bags Xj and X^. The matchings of M{i) are 
generated by combining all the pairs of matchings (M, M') G M{j) x M{k). A combination is valid if 
and only if it satisfies both the matching and cycle-free conditions. The matching condition says that a node 
cannot be matched in both M and M\ which is checked by a logical AND operation (v(M) AND w{M')). 
The cycle-free condition is checked with the union-find structures M and M': the combination is valid if 
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no node of the component of a matched node in uf (M) is a neighbor of the same component in uf(M') 
and vice versa, each test requiring 0{wf) per component. If a combination is valid, its structure M'' is 
defined by v(M'') = v(M) OR v(M'). The union-find structure is initialized from uf (M), and updated 
as the introduce bag for each matched nodes of M^ Finally, m{M") = m{M) + m{M'). As in the forget 
bag, two combinations may result in equivalent matchings, and we must compare them pairwise and choose 
the representative with the lowest number of unmatched forgotten bags. Note that the sets of forgotten 
nodes of Xj and forgotten nodes of have to be disjoint by the coherence of Definition |3] and hence no 
forgotten node can be matched twice in this setting. Furthermore, all possible combinations of matched and 
unmatched nodes are enumerated in M{j) and M{k) and hence no possible matching is overseen. 

Running time Each pair of matchings is validated and updated in time 0{wi • w'f • \og{wi)). The compar- 
ison and the choice of representative is done on-the-fly when filling the hash structure of M{i). There are at 
worst (2^^+^^)^ pairs. Thus, the complexity of the join bag dominates all other running times. Therefore, 
the complexity of the algorithm is in 0(4^^ • • \og{wi)) per bag. 

Root bag Let Xr = {x} be the root of T. M{r) consists of at most two matchings v(M) = [0] or 
v(M') = [1], where uf(M) is an empty list and \xi{M') is defined by c{x) = x. It follows that the minimum 
number of unmatched nodes for any alternating cycle-free matching of G is given by m = min{m(M) + 
l,m(M')}, and the maximum size of an alternating cycle-free matching is given by (n — m)/2 where 
n = I A^l denotes the number of nodes of G. 

Total Running Time The total time complexity of the algorithm per bag is bounded above by the running 
time of the join bag. Since there is a linear number of bags, and since for every bag Xi we have \Xi\ < 
tw{G) + 1 = It; + 1, the total time complexity of the algorithm described above is 

0(4^'+^-i(;^-logH-n). 
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Figure 8: Algorithm execution on a small bipartite cycle (top left) with its nice tree decomposition (center). 
At each bag, a set of matchings M{i) is generated according to the bag type. M{i) is represented on the 
side of each bag, with the nomenclature illustrated at the top of the figure. 



24 



F Proof of Theorem |6l 



Theorem 6. Let G be the spine of a simplicial triangulation of a 3-manifold T. Iftw{V{T)) < k, then 
tw(G) < lO/c + 9. 

Proof Let T be a tree decomposition of the dual graph, where each bag Xi contains less or equal /c + 1 
tetrahedra. We show how to construct a tree decomposition T' of the spine of T, modeled on the same 
underlying tree as T, in which each bag X[ contains less or equal 10(fc + 1) edges and triangles. 

For each bag Xi of T, we simply define the bag X'- to contain all edges and triangles of all tetrahedra in 
Xi. It remains to verify the three properties of a tree decomposition (Definition [3]). 

Node coverage It is clear that every edge or triangle in the spine belongs to some bag X[, since every edge 
or triangle is contained in some tetrahedron 5, which belongs to some bag Xi. 

Arc coverage Consider some arc in the spine. This must join a triangle t to an edge e that contains it. Let 
5 be some tetrahedron containing t\ then 5 contains both t and e, and so if Xi is a bag containing 5 then the 
corresponding bag X[ contains the chosen arc in the spine (joining t with e). 

Coherence Here we treat edges and triangles separately. 

Let t be some triangle in the simplicial complex. We must show that the bags containing t correspond 
to a connected subgraph of the underlying tree. If t is a boundary triangle, then t belongs to a unique 
tetrahedron 5, and the bags X'- that contain t correspond precisely to the bags Xi that contain 5. Since the 
tree decomposition T satisfies the connectivity property, these bags correspond to a connected subgraph of 
the underlying tree. If t is an internal triangle, then t belongs to two tetrahedra 5i and 82, and the bags X[ 
that contain t correspond to the bags Xi that contain either 5i or 82. By the connectivity property of the 
original tree, the bags containing 5i describe a connected subgraph of the tree, and so do the bags containing 
82. Furthermore, there is an arc in the dual graph from 81 to 82, and so some bag Xi contains both 81 and 82. 
Thus the union of these two connected subgraphs is another connected subgraph, and we have established 
the connectivity property for t. 

Now let e be some edge of the simplicial complex. Again, we must show that the bags containing e 
correspond to a connected subgraph of the underlying tree. This is simply an extension of the previous 
argument. Suppose that e belongs to the tetrahedra 5i, . . . , 5m (ordered cyclically around e). Then for 
each 8j, the bags Xi that contain 8j describe a connected subgraph of the underlying tree, and the bags X'- 
containing e describe the union of these subgraphs, which we need to show is again connected. This follows 
because there is a sequence of arcs in the dual graph (5i, ^2), {82^ ^3) and so on; from the tree decomposition 
T it follows that the subgraph for 81 meets the subgraph for 82, the subgraph for 82 meets the subgraph for 
(^3, and so on. Therefore the union of these subgraphs is itself connected. □ 
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